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Abstract. —We present a quantitative study of the effects of sea surface temper¬ 
atures on eelgrass productivity variables. We compared standing stock variables 
for Zostera marina for the strong El Nino event of 1986-1987, previously pub¬ 
lished by other workers for San Quintin Bay, Baja California, to our previously 
unpublished data for the “normal year” of 1992—1993. We found significant dif¬ 
ferences for most of the variables measured, which included Leaf Area Index, 
leaf length, width, dry weight and area, biomass, shoot density, and number of 
leaves per shoot. Inspection of the multivariate ENSO index (MEI) and sea surface 
temperature (SST) anomalies for each of these years showed that the differences 
could be explained by the warm SSTs associated with the ENSO event. We were 
able to explain the observed differences from a dynamic perspective by using a 
leaf-growth model forced by SSTs. We conclude that sea surface temperature 
summarizes the fundamental environmental influences on eelgrass leaf dynamics 
observed in our study site. That is, higher SSTs explain the reduction in mean 
leaf lengths and the corresponding diminution in related productivity variables. 
This study also strengthens the view that the onset of an El Nino event provides 
anticipatory evidence for the effects that a rise in global temperature is expected 
to elicit in eelgrass beds. 


Human-induced increases in atmospheric concentrations of greenhouse gases 
are expected to exacerbate global temperature change (Ramanathan 1988). The 
associated alterations could be particularly deleterious for estuaries and coral 
reefs, currently under noticeable anthropogenic stress. These relatively shallow 
environments will be affected through alterations in abiotic variables like tem¬ 
perature and sea level, wind patterns and storminess, availability of precipitation 
water and runoff, and modifications in availability of nutrients (Dyer 1985; Eman¬ 
uel 1987; Wigley and Raper 1987; Bakum 1990; Peterson et al. 1993; Watson et 
al 1996; Kennedy et al. 2002). 

Sea-surface temperature has been found to have a great influence in marine 
ecosystem dynamics (Tegner and Dayton 1987; Baumgartner et al. 1992; Beer 
and Koch, 1996; Holbrook et al. 1997; Johnson et al. 2003). In particular, warmer 
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temperatures are expected to influence the biology of organisms by reducing the 
concentration of dissolved oxygen in seawater. Diurnal and annual temperature 
oscillations also influence the biology and distribution of aquatic angiosperms 
(Holmes & Klein, 1987). Several workers have reported that growth dynamics of 
the temperate seagrass Zostera marina L. is highly correlated with sea surface 
temperature (Short an Neckles 1999; Solana-Arellano et al. 1997, Poumian-Tapia 
& Ibarra-Obando 1999; Solana-Arellano et al. 2004). According to Setchel (1929) 
eelgrass can grow only within a fixed temperature range. Other authors (Ras¬ 
mussen 1977; Phillips and Backman 1983) have also shown that temperature is 
fundamentally important in controlling the seasonal growth cycle of eelgrass. 

The distribution and abundance of seagrasses in temperate littoral waters are 
also controlled by light availability (Backman and Barilotti 1976; Denninson and 
Alberte 1985; Bulthuis and Woelkerling 1983; Orth and Moore 1988; Zimmerman 
et al 1991). In particular, light has been shown to influence the distribution (Den¬ 
nison & Alberte 1985), density (Mukai et al. 1980), flowering (Phillips & Back- 
man 1983), biomass (Mukai et al. 1980) and production (Bulthuis 1987) of Zos¬ 
tera marina. 

Dissolved oxygen, inorganic nutrients (including carbon) and water movements 
also modify photosynthesis in aquatic plants. While light, temperature and dis¬ 
solved oxygen regulate instantaneous photosynthetic rates, the availability of in¬ 
organic nutrients affects the long-term response of photosynthesis by controlling 
the levels of photosynthetic enzymes and pigments (Solana-Arellano et al. 1997). 
Eelgrass can absorb nutrients either from the roots or the leaves (McRoy et al. 
1972). Hence, modifications in upwelling activity, stratification, and tidal dynam¬ 
ics could alter the availability of dissolved nutrients and thus affect seagrass pro¬ 
ductivity. In our study area (San Quintin Bay, Baja California, one of the last 
remaining ecologically functional estuaries in the California Bight), the onset of 
an ENSO event alters sea surface levels and temperatures, as well as tidal dy¬ 
namics and nutrient availability (Alvarez-Borrego 2004). Hence, it is reasonable 
to expect that ENSO-like conditions would influence eelgrass dynamics in the 
bay. Sea surface temperature provides a reasonable paradigm for an ENSO-like 
alteration. Since thermal expansion of the water column influences sea level var¬ 
iation (Wigley and Rapper 1987) sea level can be assumed to increase, thus mod¬ 
ifying light availability and changing tidal dynamics and circulation (Short and 
Neckles 1999). Also, thermally induced stratification would imply a deeper ther- 
mocline and nutricline, and modify upwelling transport of dissolved nutrients 
(Petterson et al. 1993). Hence we would expect sea surface temperature to be a 
variable that triggers changes in the dynamics of the principal environmental 
influences on eelgrass in our study site. 

To provide a basis for the assumption of significant causal linkage between an 
ENSO warming event and changes in the growth dynamics of eelgrass we dem¬ 
onstrate here that sea surface temperature accounts in a fundamental way for the 
observed variability of standing stock variables by comparing two whole-year- 
cycle data sets of Z. marina L. in the San Quintin Bay estuary. The first data set, 
for the ENSO event of 1987-1988, was previously published by Ibarra-Obando 
et al. (1997). The second data set was collected by our research group from 
November 1992 through November 1993. Both data sets consist of means for leaf 
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length, width, dry weight and area, biomass and density per square meter, Leaf 
Area Index (LAI), and number of leaves per shoot. 

Due to the architecture of Z. marina , above-ground productivity is determined 
by the dynamics of leaf elongation. To show that temperature controls productiv¬ 
ity, we demonstrate a causal relationship between temperature and leaf length by 
using a leaf-growth model forced by sea surface temperature. We produced a 
consistent representation of leaf-length dynamics for both data sets and corrobo¬ 
rated that sea surface temperature can be considered a fundamental environmental 
forcing agent for productivity variables of eelgrass at our study site. An alteration 
of the normal variation range for this variable like that registered during the onset 
of an ENSO event provides anticipatory evidence for the effects that an equivalent 
change in global temperature would elicit on eelgrass beds at our study site. This 
agrees with the view on this causal relationship for seasgrasses discussed in Short 
and Neckles (1999). 


Study Site 

San Quintin Bay is a Y-shaped shallow coastal lagoon on the Pacific side of 
the Baja California Peninsula, Mexico (30 30N-116°10'W), which has a total area 
of 42 km 2 . Evaporation exceeds runoff plus precipitation. In the Bay’s waters, 
tides are mixed and predominantly semidiurnal and are considered the main cause 
of temperature variability (Alvarez-Borrego 2004). Upwelling events bring nutri¬ 
ent-rich waters near the mouth of the bay (Dawson 1951) and tidal currents prop¬ 
agate those waters throughout the bay. To some extent, this accounts for the rather 
high nutrient content of its waters (Barnard 1962; Alvarez-Borrego and Alvarez- 
Borrego 1982). Organic matter is trapped and materials are re-mineralized, re¬ 
leasing nutrients back to the water column. These remineralization processes 
(Smith et al. 1991), along with turbulence induced by tidal currents and waves, 
control nutrient concentrations in the lagoon (Alvarez-Borrego 2004). 

Sampling and Exploratory Analysis 

The sampling site is a mudflat located in the west arm of the bay. We collected 
shoots monthly from sixteen 20 X 20 cm quadrats selected at random, from 
November 1992 through November 1993. Each sample was placed in individually 
labeled plastic bags and cooled until processed. Each shoot was cleaned with 
distilled water. We determined the number of shoots and leaves per unit area, the 
length (mm), width (mm) dry weight (g), and area (mm 2 ) of each leaf. A total of 
10,000 complete leaves were collected through the entire year cycle. Biomass per 
unit area and L.A.I estimations are presented in g m~ 2 and m 2 m -2 respectively. 
For both data sets that were analyzed, underwater radiation at 1 meter under sea 
level was measured (Em 2 ) by using a PAR and direct beam quantum radiometer- 
photometer (Li-Cor, Inc.) using an integration time of 1000 s. The integrated 
photon flux rate in a time interval centered at noon gives a reasonable estimate 
of the maximum radiation in the day. Daily variation can then be approximated 
by the Monteith sine law (Monteith 1965). In situ temperatures for both data sets 
were also obtained (figure 5 b). Principal Component Analysis (PCA) was de¬ 
veloped for all variables in both our data and that of Ibarra-Obando et al. (1997) 
to assess the importance of abiotic variables. Nutrient data are not available for 
1992-1993 so we used upwelling index values as a proxy in accordance with 
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Alvarez-Borrego (2004) that reports upwelling as the main source of dissolved 
nutrient availability in the site. 

Data in MEI values SST anomalies and Upwelling index were obtained from 
NOAA (http://www.cdc.noaa.gov.people.klauss.wolter/MEI/table.html, http://ferret. 
pmel.noaa.gov/NVODS/servlets/datasets and ftp://orpheus.pfge.noaa.gov/outgoing/ 
upwell/monthly/upindex.mon). Data on sea level values at San Quintin Bay were 
provided by the Laboratorio de Nivel del mar, Departamento de Oceanografia 
Fisica CICESE (http://nivelmar.cicese.mx/). 

Theoretical Methods 

Following Nadezhda et al. (2001) we considered that water temperature, light 
radiation and nutrients are the fundamental environmental influences on eelgrass 
growth. In accordance with Alvarez-Borrego (2004) we will consider that the 
availability of dissolved nutrients is controlled by upwelling activity at our study 
site. We designate sea surface temperature by T(t) at a time t. Similarly cf>(0 means 
underwater light radiation and U(t) are values of the upwelling index. 

Let l(t) denote eelgrass leaf length. Then, leaf length rate dynamics will be 
assumed to follow the growth model, 

“ = *«[/- - /(OL (1) 

clt 

where l x is the maximum possible value that leaf length can be attain. This con¬ 
stant is also known as the asymptotic upper bound for leaf growth. For leaf length 
values that are very close to L, vanishing leaf length rates are expected (Batschelet 
1974). The scaling factor k (t) summarizes the fundamental environmental influ¬ 
ences on leaf dynamics. It can be formally represented through the expression 

k(t) = f(T(t\ 4>(f), um (2) 

where/(0, <J>(0* ^(0 * s assumed to be continuous and differentiable in all of its 
arguments. Increasing water temperatures are known to affect eelgrass metabolism 
and the maintenance of a positive carbon balance (Zimmerman et al. 1989). Op¬ 
timum eelgrass photosynthesis is attained at temperatures below a seasonal max¬ 
imum (Biebl and McRoy 1971). Temperatures above the growth optimum to near 
the upper limit of thermal tolerance have resulted in reduced eelgrass productivity 
(Thayer et al. 1984; Moore et al. 1996). This conceptual framework, along with 
the assumption that a change in T(t) triggers a response in the remaining variables 
of equation (2) is consistent with the assumption that k(t) can be represented in 
terms of mean sea surface temperature T(t). The forcing factor k(t) can be assumed 
to have the polynomial form 

N 

k(t ) = 2 a n T"(t). (3) 

n = 1 


Integration of equation (1) from t t to t i+l led us to the equation, 

ly, ~ l(t;+ 1) 


exp 


r 


k(t) dt 


l - m 

After a few algebraic steps equation (4) provides the regression expression. 


( 4 ) 
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Fig. 1. Mean monthly values for Zostera marina standing stock variables and leaf area index 
(LAI) measured in San Quintin Bay during the “normal” year of 1992-1993. 


l(t i+ ,) = L~ IL ~ (5) 

where by virtue of equation (3) G(/„ ) has the form: 

G(?„ t M ) = -Ij cij J * T'(t) dt (6) 

In order to corroborate empirically the assumption of equation (3) we fitted 
equation (5) to available data on leaf length and sea surface temperature. Deter¬ 
mination coefficients, graphs of predicted vs observed values and residual analyses 
are presented in the results section. The mean square regression method and the 
statistical package STATISTICA were used in these fittings. 

Results 

Maximum leaf width was 4.6 mm and occurred in November 1992, and max¬ 
imum leaf length of 290 mm was reached in August 1993 (Fig 1). Leaf-widths 
close to maximum values also occurred in September and November 1993. 
Changes in mean leaf length and width followed the same variation pattern; the 
maximum for both measurements was attained by late summer (Fig. la). Similar 
results were obtained for mean leaf area and mean leaf biomass; both measure¬ 
ments reached their maximum in August, (Fig. lb), these were 1543 m m 2 and 
0 .034 g respectively. Biomass per unit area reached a maximum of 57.02 g m 2 
in August, but another peak of 54.04 g m 2 occurred in May (Fig. lc). This was 
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Fig. 2. Comparison of mean values for variables in our study with those of the 1986-1987 El 
Nino year reported by Ibarra-Obando et al. (1997) data at the same depth. Dashed lines corresponds 
to our study and continuous lines to Ibarra-Obando et al. (op. cit) 

likely due to the period of reproduction, as a large number of reproductive shoots 
was observed at this time. Mean L.A.I had a maximum of 3.3 and was also 
attained in May. The inverse relationship between mean numbers of shoots and 
leaves is shown in Fig. Id, such that minimum numbers of shoots (less than 400 
shoots per square meter), occur when the number of leaves reaches its maximum 
of 6 leaves per shoot, and maximum number of shoots (1390 shoots m 2 ) is reached 
almost at the same time the average number of leaves per shoot is declining to 
its minimum (2 leaves per shoot). This relationship could be explained by the 
autoecology of seagrasses (Kentula and Mclntire, 1986), since light decreases 
during winter, the plant uses its energy to produce fewer shoots with a large 
number of leaves so they have more access to light. 

For the El Nino year of 1986-1987, Poumian-Tapia & Ibarra-Obando (1999) 
reported smaller values for mean length, biomass, number of shoots and number 
of leaves per shoot than the ones found in our study for the same area (Fig 2.). 
We fitted a general linear model for upwelling index values in terms of dissolved 
nitrates and phosphates and their cross products, finding a determination coeffi¬ 
cient of R 2 = 0.88. All estimated parameters had p-vcilues < 0.05 ( 0.0090 for 
nitrates, 0.009 for phosphates and 0.03 for the cross product of nitrates and phos¬ 
phates), meaning that the 3 terms contribute to the multiple correlation coefficient. 
This result also provides a quantitative basis for the linkage between dissolved 
nutrient concentration and upwelling index values. 
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A Principal Components Analysis (PCA) for the 1987 data extracted a first 
component defined mainly by temperature, with a correlation coefficient of R = 
0.95 and irradiance with a correlation coefficient of R = -0.74. This component 
explained 39% of the variability of biotic variables. Second and a third compo¬ 
nents were represented by nitrates and phosphates in that order, with coefficients 
of correlation R = 0.91 and R = -0.66 which explained 29% and 20% of the 
variability, respectively. 

For the 1992-1993 data we used upwelling index values as a proxy for dis¬ 
solved nutrient availability. In that PCA, we found a first principal abiotic factor 
defined by temperature and irradiance with correlation coefficient of R = 0.85 for 
temperature and —0.68 for irradiance. This factor explained 41.4% of the vari¬ 
ability of biotic variables. A second component was mainly defined by upwelling 
with R = 0.74 and this factor explains 33.2% of the variability of all variables. 

PCA results for both years, which showed a dominant forcing of SST on eel- 
grass productivity variables, were corroborated by the fittings of equation (5). 
This fitting produced a high coefficient of determination (R 2 = 0.94) for the 1987 
data and an R 2 = 0.97 for the corresponding 1992-1993 leaf length data. Analysis 
of residuals showed a good correspondence between eelgrass leaf length and sea 
surface temperature for both data sets. Figure 3 shows the corresponding predicted 
versus observed values associated for each fit. This demonstrates that the observed 
differences can be explained by higher temperatures. In other words, sea surface 
temperature is a variable that summarizes the fundamental environmental influ¬ 
ences that determine leaf length dynamics, and therefore of other variables such 
as biomass and L.A.I for eelgrass in our study site. 

Discussion 

This work presents the most recent data for important variables used in eelgrass 
production studies in San Quintin Bay (fig. 1). The nearest related sampling in 
the area occurred in 1987-1988, as reported by Ibarra-Obando et al. (1997). They 
reported smaller values for mean leaf length and also for most of the productivity 
related variables than our present results (fig. 2). An analysis of the causal linkage 
between sea surface temperature, underwater light radiation, dissolved nutrient 
availability and eelgrass productivity variables leads us to conclude that a higher 
temperature is the most important factor determining the observed differences. 
Our findings sustain the paradigm that a warm temperature stress can induce a 
deleterious effect in Z. marina , (Rasmusen, 1977, Penhale, 1977; Wetzel and 
Penhale, 1983; Evans et. al. 1986; Johnson et al. 2003). 

Primary and secondary production in the California Current (CC) are high when 
the transport of the cool, low-salinity water from the north is strong and low when 
this transport is weak (McGowan 1983). Periods of warming and low productivity 
in the CC tend to coincide with El Nino events in the eastern tropical pacific. 
ENSO events induce warming of SSTs in San Quintin Bay (Silva-Cota and Al- 
varez-Borrego, 1988). The present analysis shows that slower growth rates in Z. 
marina reported by Ibarra-Obando et al. (1997) could be a consequence of the 
ENSO event that took place from June 1986 to January 1988, just before and 
during the year where their data were taken Comparison of ME1 and SST values 
demonstrates clear differences between the 1987 ENSO and the “normal” con¬ 
ditions of 1993 (Fig 4.). 
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Fig. 3. Observed versus predicted values for the fitting of equation (5) to mean length data, a) 
Fitting corresponding to the Ibarra-Obando et al. (1997) data, b) Fitting corresponding to our data. 


Thermal expansion of sea water results in higher sea level values (Wigley and 
Rapper 1987), which can reduce light reaching seagrass beds at depth and thus 
limit photosynthesis (Short and Neckles (1999)). For Z. marina, the effects of 
decreased light are reduction in shoot density, leaf width, number of leaves per 
shoot, growth rate (Short et al. 1993) and above-ground plant productivity (Back- 
man and Barilotti 1976). Comparison of sea-level anomalies for San Quintin Bay, 
shows largely negative anomalies during the nQn-ENSO conditions of 1985 and 
positive anomalies during the 1987 ENSO (fig. 5). Thus we can infer that sea 
level anomalies were largely negative in 1993. There were also significant dif- 
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Month 

Fig. 4. Above) Multivariate Enso Index values for 1987 and 1993. Below) Corresponding Sea 
Surface Temperature anomalies. 


ferences (p = 0.03) in mean light radiation between our 1993 data and those of 
Ibarra-Obando (1997). This occurred for both the incident light flux and under¬ 
water light availability (Fig 6). Temperature stress and reduced availability of 
light likely contributed to lower productivity of eelgrasss in 1987. 

Upwelling of nutrient-rich waters during non-ENSO conditions maintains high 
nutrient levels in the water column. But thermal stratification can modify the depth 
of both the thermocline and the nutricline and so influence nutrient transport by 
upwelling. Upwelling index values are in a good correspondence with dissolved 
nutrient availability, corroborating reports on the role of an ENSO event in the 
reduction of dissolved nutrient availability in our study site (Silva-Cota and Al- 
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Fig. 5. a) Sea level anomalies for San Quintin Bay for the ENSO year of 1987 and 1985 a year 
before the onset of the referred ENSO event. Higher sea level values for the ENSO event year area 
observed, b) Mean monthly temperatures for Ibarra-Obando et al. (1997) and present study data sets 
at San Quintin Bay. 


varez-Borrego, 1988). This strongly suggests an unfavorable influence for eelgrass 
productivity, in addition to temperature and light availability effects. 

We conclude that for the 1987 data the conjunction of a high temperature stress, 
reductions in light radiation and availability of dissolved nutrients could explain 
the diminution in eelgrass productivity variables during the 18-month ENSO 
event. Ibarra-Obando et al. (1997) found a turnover time of about 38 days for Z. 
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Fig. 6. Above) Incident Light Radiation in San Quintin Bay during the ENSO event year of 1987 
and 1993 which was consider as normal. Below) Underwater light luminosity corresponding to the 
same years. We can notice a grater light availability for the 1993 sampling period. 


marina L., in San Quintin Bay. Reduced turnover times could have shortened 
responses to environmental stress for these seagrass populations in such a way 
that the associated effects were observable during the reported year—cycle study. 
In any event, both the performed PCA and the fitting of equation (5) show that 
environmental influences which controlled the observed dynamics can be sum¬ 
marized in a reasonable way by sea surface temperature variability. 

The present study provides the first quantitative framework for corroborating 
the assumption that an increased temperature, as observed on the onset of an El 
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Nino event, will elicit reduced growth rates in Z. marina (Fig. 3). ENSO driven 
alterations are limited to 12-18 months between periods of normal conditions that 
range from 2-7 years. A continued ENSO-like anomaly could result in long-term 
unfavorable temperature conditions, triggering changes in light radiation and in 
the availability of dissolved nutrients. Besides the abatement of productivity var¬ 
iables, this could also disrupt adaptive response in eelgrass and promote its re¬ 
placement by subdominant species as reported by Johnson et al. 2003. 

Global climate change is expected to result in modifications of SST, sea level, 
C0 2 concentrations and UV radiation reaching the Earth’s surface (Houghton and 
Woodwell, 1989; Mitchell 1989; Schneider 1994; Kerr 1992; Watson et al. 1996; 
Titus 1990; Smith et al. 1992; Schellnhuber et al. 2006), and to induce severe 
effects on both terrestrial and aquatic plants (Watson et al. 1996). In particular, 
global change is likely to exacerbate the deleterious effects of human activities 
on estuaries and costal lagoons (Kennedy et al. 2002). The effects of anthropo¬ 
genic activities on estuaries are well documented in the literature (Roblee et al., 
1991; Walker and McComb 1992; Short and Willie-Echeverria 1996). On the 
other hand the mechanism underlying climatic change effects are poorly under¬ 
stood. The expected changes linked to global warming on submerged plant species 
like sea grass beds and salt marshes are of a great concern due to the important 
ecological services they provide and the current extent of their loss. 

The present study indicates that the onset of an El Nino event alters the San 
Quintin Bay environment with changes in sea-surface temperature, level and un¬ 
derwater light radiation. Similar changes for these variables are expected to occur 
in a global change scenario. Even though ENSO events have also driven changes 
in dissolved nutrient availability at the site, it is more difficult to predict how 
global change will affect this factor. Upwelling caused by northerly winds and 
the resulting offshore Eakman transport is a dominant oceanographic process in 
spring and summer along the entire California and Baja California Coast (Sver¬ 
drup et al. 1942). Observations show that wind-driven upwelling along the Cal¬ 
ifornia Current has increased over the past 30 years (Snyder et al. 1993). This 
increase is assumed to be elicited by temperature gradients induced by increasing 
greenhouse gas forcing, but such association has been speculative. In the event 
of a positive association of global climatic change and intensification of upwell¬ 
ing, it might be reasonable to expect that the observed non-limiting nutrient en¬ 
vironment within San Quintin Bay could be maintained. But increased sea tem¬ 
perature in the waters adjacent to the mouth of the bay and a coupled stratification 
could promote a deepening of the thermocline and nutricline, rendering upwelled 
waters deficient in nutrients, regardless of how intense the transport might be 
(Petterson et al. 1993). This coincides with reports on a reduced dissolved nutrient 
availability in the waters adjacent to the mouth of San Quintin Bay during an 
ENSO event (Silva-Cota and Alvarez-Borrego 1988). 

Our findings strentengthen the view that the onset of an El Nino event provides 
valuable anticipatory scientific evidence on the effects on eelgrass productivity 
variables as a consequence of a global climatic change (Schneider 1994), partic¬ 
ularly those associated with sea surface temperature, and underwater light radia¬ 
tion. For dissolved nutrients, the onset of an El Nino event at our study site could 
reduce dissolved nutrient availability probably linked to a deficient upwelling 
transport. 
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